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Abstract 

We present a connection between the ground state of the Ising 
Model and the dimer problem. We find the generating function for 
dimers as the appropriate limit of the free energy per spin for the Ising 
Model. 

1. Dimers. Dimers are objects connecting two neighboring sites of the 
lattice. We consider a lattice which can be covered with dimers in the way 
that each site of the lattice belongs to exactly one dimer. How to calculate the 
number of such coverings for a given lattice was shown in 1961 by Kastelayn, 
Fisher and Temporley in 1961 [1-3]. The connection between the Ising model 
on the plane lattice and the dimer problem is well known. In 1961, Fisher 
[4] gave a general rule connecting the partition function for the Ising model 
with a generating function for configurations of dimers on an appropriately 
constructed lattice. For the Ising model on the honeycomb, for instance, the 
dimers on 3-12 lattice are considered. Moreover, the Ising model — as well 
as the model of free fermions — are connected with the dimer problem on a 
nonplane lattice [5,6]. 

2. From the Ising model to the dimer problem. In this paper, we 
shall show that the doubled number of dimer coverings for a given lattice is 
equal to the degeneration of the ground state of the Ising model on a lattice 
dual to the one covered with dimers. Next, we will use this fact to find the 
generating function for the dimers on the chessboard lattice. 

Let us begin with the description of our construction. Let L be any 
plannar lattice for which the dual lattice can be covered with dimers. Let 
us cover the lattice with dimers in one of the possible ways. 



1 



To every edge crossed by a dimer we assign the positive constant of inter- 
action E, and to the remaining edges — the negative constant of interaction 
—E. 

The energy of the interaction of two neighboring spins is defined by the 
product of these spins and the constant of interaction. The degeneration of 
the ground state equals the number of spin configurations with the lowest 
energy. 

The degeneration of the ground state of our Ising model is equal to half 
of the number of the dimer coverings of the lattice Lj. 

To prove the above statement, let us note that we get the smallest energy 
when exactly one edge of each elementary polygon of the lattice L has the 
positive energy E. 

Each configuration of this type corresponds to one dimer covering of the 
lattice Ld- Conversely, each dimer covering corresponds to two spin config- 
urations with the lowest energy. To see that, we fix one spin and systemat- 
ically define the neighboring ones: if the edge is positive — or negative but 
crossed — the next spin has the same sign; otherwise the next spin has the 
opposite sign. 

Obviously, we can fix the first spin in two ways. Hence, the number of 
dimer coverings reads 

Z d = - lim Z(T)e ME / kT , (1) 

where Z is the statistical sum of the Ising model on the lattice L, ME is a 
minimal energy of the system, and M is the number of edges of the lattice 
L minus number of dimer s. 

3. Example. Finally, we present an application of the procedure described 
above. Let us consider the chessboard Ising model with constants of interac- 
tion E x , E 2 , E 3 <0,E 4 >0 (see Fig.l). 

We change the signs of E 2 and E 4 in every other column so that the edges 
crossed by dimers of the standard configuration have positive interaction 
constant (Fig.l). This transformation does not change the statistical sum. 

Let Ei = -E - LxkT, E 2 — —E — L 2 kT, E 3 = —E — L 3 kT, E 4 = 
E + L 4 kT. 

Then 

lim Z(E 1 ,E 2 ,E 3 ,E 4 )e ME / kT = Z d (x 1 ,x 2 ,x 3 ,x 4 ), (2) 
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where Z(xi,x 2 ,x 3 ,X4) is the generating function for the dimer problem on 
the chessboard lattice in the variables Xi = exp(2Lj — LI — L2 — L3 — LA), 
M = 3N/4, and N is the number of spins. 

In 1951 Utiama obtained the formula for the free energy per spin for the 
Ising model on the chessboard lattice [6]: 

f(T) = -kT hm = 

= -kT\n2-— / hi -(1 + C 1 C 2 C 3 C A + S 1 S 2 S 3 S A + (3) 
Lon z Jo Jo 2 

+ (SiS 2 + S3S4) COS0 - (S1S4 + S 2 S 3 ) COS0 + 

+S 2 S 4 cos(0 + 6) + S1S2 cos(0 - 6))d6d(i>, 

where d = cosh(-2^/A;T), and $ = sinh(-2£;/A;T). 

While calculating the thermodynamic limit, we suppose that the lattice 
grows in the same way in two directions. 

Using this result we find 

*{ Xl ,x 2 ,x 3 ,x 4 )- Jim - _- r hm— - (4) 



j /■27T p2tt 

— / / \n(xj + x\ + x\ + xl+ 
bir z Jo Jo 



16?r 2 

+2(a; 3 a;4 — X\X 2 ) cos0 + 2(x 2 x 3 — X1X4) cos9— 

—X\X 3 cos(0 + 6) + a^a^ cos(0 — 9))d9d(f>. 

In this way it was possible to find the generating function for the dimer 
problem without any refernce to the general theory. Similar considerations 
may be repeated for the triangular and honeycomb lattices. 
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Fig.l. Lattices: L (continues lines) and L d (dotted lines). 
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